1/1 


UNCLASSIFIED 


ELECTROMAGNETIC DESCRIPTION OF THE THO STREAM 
INSTABILITV IN A RELATIVISTIC ELECTRON BEAM(U> NAVAL 
SURFACE UEPPONS CENTER SILVER SPRING MD H S UHM MAV 84 
NSHC/MP-84-236 F/G 28/8 NL 



































AD-A179 864 


.» J* V ~-P "J« ” V - TT" 


r vn -* v « ITT^. ’ V^"WV, l"rrr r* » V » V^T:v r ;y_ , y. ,, vry. " V l V 



ft 


IfflLflLLCOPJf 

ELECTROMAGNETIC DESCRIPTION OF THE TWO 

STREAM INSTABILITY IN A RELATIVISTIC 
ELECTRON BEAM 


I r 

ft.. 


HAN S. UHM 



NSWC MP 84-236 

PLASMA PHYSICS PUBLICATION NO. 84-4 
MAY 1984 




NAVAL SURFACE WEAPONS CENTER 

White Oak, Silver Spring, Maryland 20910 


V ^ 




ELECTROMAGNETIC DESCRIPTION OF THE TWO STREAM 
INSTABILITY IN A RELATIVISTIC ELECTRON BEAM 



Han S. Uhm 

Naval Surface Weapons Center 
White Oak, Silver Spring, Maryland 20910 


The influence of electromagnetic effects on the two stream instability in 
a relativistic electron beam propagating through a collisionless plasma 
channel is examined, and closed algebraic dispersion relation for the complex 
eigenfrequency o> is obtained. It is shown that electromagnetic effects can 
have a strong stabilizing influence on the two stream instability in a 
relativistic electron beam, drastically enhancing the critical beam current 

for instability. 
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One of the most basic instabilities that characterize a relativistic 
electron beam propagating through a collisionless plasma is the two stream 
instability, 1 " 3 which resulted from the relative drift motion between the oeam 
electrons and the background plasma particles. Although the two stream 
instability is relatively familiar in the plasma physics community, most of 
the previous studies 1 ’ 2 on this instability have been limited to one 
dimensional calculations. In recent literature, 3 Bogdankevich and Rukhadze 
investigated the two stream instability of a relativistic electron beam, 
including the finite radial geometry effects on stability behavior. Thus, 
they were able to determine the limiting beam current due to the two stream 
instability. However, their calculation has been based on the electrostatic 
approximation. Although this is a reasonable approximation for a mildly 
relativistic electron beam, a significant modification to the stability 
behavior is expected for a ultrarelativistic beam; where the electromagnetic 
effect often plays an important role. In this letter, the influence of 
electromagnetic effects on the two stream instability is investigated for a 
relativistic electron beam propagating through a collisionless plasma channel. 

The analysis in this paper is carried out within the framework of a 
macroscopic cold fluid model assuming that the beam-plasma fluids are immersed 
in a uniform axial magnetic field. For purposes of analytic simplification, 
the stability analysis is specialized to the case of a sharp-boundary 
equilibrium in which the beam and plasma channel have rectangular density 
profiles, i.e.. 
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0 < r < R b , 
















where the subscript j = b, i, e denote beam electrons, plasma ions and 
electrons, respectively, r = R c is the radial location of a grounded 
conducting wall. Since the unstable mechanism of the two stream instability 
is mostly due to the fluctuations of the axial electric field, the stability 
analysis in this article is restricted to perturbations of the transverse 
magnetic mode polarization. In addition, it is further assumed that the 
perturbations are axisymmetric (3/30 = 0). This is a reasonable assumption 
because the axisymmetric perturbation is the most unstable mode for this kind 
of instability. Combining all of these assumptions and restrictions, it is 
straightforward to show that the eigenvalue equation of the two stream 
instability is expressed as 
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where 6E (r) is the perturbed axial electric field, 5p(r) = l e.6n. and 
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<5J (r) = l eilrj <5n. + nVsv- ) are the perturbed charge and the axial 

4 : J J J J ^ 

J 0 
component of the current densities, respectively, e. is the charge, n. is the 
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equilibrium density, v9 is the mean equilibrium axial velocity, 6n. is the 
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perturbed density and 6V. is the perturbed axial velocity, of species j, w is 
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the complex eigenfrequency and k is the axial wavenumber. 


Consistant with polarizations of the two stream instability, the 
transverse component of the perturbed electromagnetic field is neglected on 
the present analysis. This approximation tremendously simplifies the 
linearized fluid calculation of the equation of motion and continuity equation 
for the beam and plasma fluid element. Thus, after a straightforward 
calculation, the perturbed axial velocity is expressed as 
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where iiij is the rest mass of particles of species j. Similarly, the perturDed 
density is given by 
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More detail calculation including the transverse component of the 
electromagnetic fields are currently under investigation by the author and 
will be presented elsewhere. Note from Eqs. (3) and (4) that both the 

A A 

perturbed density 6n. and axial velocity 6V- profiles can eventually be 
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expressed in terms of the perturbed axial electric field 6E Z . Substituting 
Eqs. (3 } and (4) into Eq. (2) and carrying out a straightforward calculation, 
the eigenvalue equation in Eq. (2) is expressed as 
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where ur.(r) = 4nnV(r)e^/m. is the nonrelativistic plasma frequency-squared of 
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particles of species j. Note that urjr) = w . = const., for 0 < r < R b and 


zero, otherwise. 


For the present purposes, the eigenvalue equation (5) is solved for long 
wavelength perturbations characterized by 
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Within the context of Eq. (6), the physically acceptable solution to Eq. (5) 
can be approximately given by 


1, 0 < r < R. , 

6E(r) -a D (7) 

' -tn{r/R c )/tn(R c /R b ), R b < r < R £ , 

where a is a constant and R c > > R b is assumed in order to satisfy Eq. (6) 
self-consistently. Multiplying Eq. (5) by r6E z (r) and integrating over r from 
r = 0 to r = R c , it is straightforward to obtain the dispersion relation 
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where use has been made of Eq. (7). In Eq. (8), or = 4ne n^/m is the 
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nonrelativistic plasma frequency-squared for beam electrons, f g = n e / n b . and 
the ion contributions are neglected. In obtaining Eq. (8), it is assumed that 
the beam electrons have the same axial velocity B b c which is related to the 

o p "1 

mass ratio y b by = (1 - B b ) . Axial velocity of the plasma electrons is 
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assumed to be zero. Within the context of Eq. (6), the term k - ur/c in the 
right-hand side of Eq. (8) can be neglected. However, I keep this term for 
the clarity in the subsequent stability discussion, without affecting the 
final results in this letter. 

Two points are noteworthy from Eq. (8). First, it is obvious that the 
term 2/R b *n(R c /R b ) in Eq. (8) represents the influence of finite radial 
geometry. Somehow, neglecting this term recovers the familiar one dimensional 
dispersion relation for the two stream instability. Second, the influence of 
electromagnetic effects are incorporated into the terms proportional 
to u /c in Eq. (8). Neglecting these terms gives the dispersion relation 
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for the electrostatic approximation. Paralleling the standard method 


developed in the previous literature, ’ the critical beam current for 


i istability is analytically found from Eq. (9) and is given by 
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for the electrostatic approximation, where I A * mc J /e - 17000 ampere is the 


Alfven critical current. The electron beam with current below the critical 


current is stable whereas the beam with current above the critical current is 


unstable. Note from Eq. (10) that the critical current is inversely 


proportional to the plasma electron density(f„). 


In order to demonstrate the influence of electromagnetic effects on 


stability behavior, Eq. (8) is numerically solved for complex eigenfrequency 
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u = o> r + iui.j, determining the stability boundaries in (w,y b ) parameter space 


for specified values of f e . The critical current I cm of the electromagnetic 


calculation is obtained from these stability boundary curves. For convenience 


in the subsequent analysis, we define the electromagnetic current enhancement £ 
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which is a function of y b and f e in general. Shown in Fig. 1 is the 
electromagnetic current enhancement £ versus the beam energy y b obtained from 
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Eqs. (8) and (11) for R c /R b = 7.4, and f g = 1 and f g = 10. Obviously from 
Fig. 1, the critical current for the electrostatic approximation is valid only 
for a mildly relativistic electron beam. The electromagnetic current 
enhancement £ increases drastically with increasing values of energy Y b and 
plasma density f g , thereby exhibiting a substantial enhancement of the 
critical current by the electromagnetic effects for a relativistic beam. 


It is instructive to rewrite Eq. (8) by 
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which is plotted in Fiy. 2. When the value of the local minimum point labeled 

P is less than unity, Eq. (12) has six real roots and the system is stable. 

It is noted from Eq. (12) that when $ b approaches to unity, the strength of 
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the destabilizing term proportional to u> p /Y b reduces drastically in comparison 
with the stabilizing term proportional to 2c £ /R b tn(R c /R b ). In Fig. 2, the 
influence of vertical line w = kc drastically pull down the local minimum 
point P as the line w = ke b c approaches to u = kc, eventually providing strong 
stabilizing effects. The stabilizing influence of electromagnetic effects 

A 

originates from the perturbed axial current density 60 b (r) in Eq. (2). 


This research was supported in part by the Independent Research Fund at 


the Naval Surface Weapons Center and in part by the Defense Advanced Research 
Project Agency. 


























REFERENCES 


Krall, N. A., and Trivelpiece, A. W., Principles of Plasma Physics , (New 
York: McGraw-Hill, 1973), Chapter 9. 

Davidson, R. C., Theory of Nonneutral Plasmas , (Mass.: Benjamin Reading, 
1974), Chapter 2. 

Bogdankevich, B. t., and Rukhadze, A. A., Soviet Physics USPEKHI, 14 , 163 























FIGURE CAPTIONS 






































































